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The quantum phase transition in the ground state of the extended spin S = 1/2 XY model has 
been studied in detail. Using the exact solution of the model the low temperature thermodynamics, 
as well as the ground state phase diagram of the model in the presence of applied uniform and/or 
staggered magnetic field are discussed. 

PACS numbers: 71.27.+a- Strongly correlated spin systems; 75.10.Jm - Quantized spin models. 



I. INTRODUCTION 



There is a considerable interest in models of strongly 
correlated electron and spin systems showing a Quan- 
tum Phase Transition (QPT) (see the Ref. 1 ) . In the 
case of One-Dimensional electron or spin systems QPT 
related to the dynamical generation of a charge or spin 
gap is often connected with the change in the topology 
of the Fermi surface, in particular with the doubling of 
the number of Fermi points 2-5 . 

The one-dimensional spin 5=1/2 XY chain 



Hxy = - J 



I QX OX 



+1 



qv qv \ 



(1) 



is the simplest exactly solvable strongly correlated spin 
model 6,T . Its exact solution is expressed in terms of the 
Fermi gas of spinless particles (spinless Fermions (SF)). 
The free Fermi gas principle for construction of the sys- 
tem eigenstates and eigenvalues provides the straightfor- 
ward and easy way to obtain the exact expressions for 
correlation functions and thermodynamic quantities (see 
the Ref. 8 and references therein). 

More than thirty years ago M. Suzuki proposed the 
whole class of generalized XY models with multi-spin 
interaction, allowing the exact solution in terms of the 
Fermi gas of spinless Fermions 9 . In the fermionic 
representation, the multi-spin coupling shows itself only 
through the form of the single particle spectrum. In this 
paper we consider the simplest of the proposed general- 
ized XY models - the extended XY model with three 
spin coupling. The Hamiltonian of the model is given by 



- J 



( cx cx 
V a n a n+2 



(2) 



and describes the spin system determined on the zig-zag 
chain (see Fig.l), with the transverse exchange between 
the spins on the nearest-neighbor sites J, and trans- 
verse exchange between the spins on the next-nearest- 
neighbor sites J*, the latter depends on "z" orientation 
of the spin being between the next-nearest-neighbors. 




FIG. 1. Schematic representation of the structure of the 
extended XY model. 

Despite its rather formal and possibly even nonrealis- 
tic nature the model (2) attracted our attention owing to 
its several important advantages: 

• it is an exactly solvable model; 

• the exact solution is expressed in terms of Fermi 
gas of SF. 

Moreover, as it is shown in this paper, the model (2) is 
characterized by the rich ground state phase diagram. In 
particular 

• with the increase of three spin coupling, at J* = 2 J 
the system experiences QPT from the Spin Liquid 
I phase into the Spin Liquid II phase; 

• in the SF representation this QPT is associated 
with the doubling of Fermi points; 

• at the transition point the magnetic and low- 
temperature thermodynamic properties of the sys- 
tem show a well pronounced anomalous behavior; 

• in the case of applied uniform and/or staggered 
magnetic field the system is characterized by the 
rich ground state phase diagram which contains 
ferromagnetic, antiferromagnetic and three differ- 
ent spin liquid phases. 

The paper is organized as follows. In the forthcoming 
Section we discuss the spinless Fermion representation of 
the model. In the Section III the ground state proper- 
ties are discussed. In the Section IV the low-temperature 
thermodynamics of the system is considered. In the Sec- 
tion V the ground state phase diagram of the model is 
studied in the presence of magnetic field. Finally, the Sec- 
tion VI contains the discussion and concluding remarks. 
The paper contains also two appendices. In Appendix 
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A we present the table with exact expressions for the so 
called emptiness formation probability (EFP) P(n) for 
n = 1, 20 and for different values of the parameter J*. 
In Appendix B we present the expression for the EFP, 
which fits our exact data for a> a c . 



II. SPINLESS FERMION REPRESENTATION 

The Hamiltonian (2) can be diagonalized by means of 
the Jordan- Wigncr transformation 10 




-1.5 



n-l 



m=l 
n-l 



$n — &n — c « (l ^ C ln C m) 



(3) 



m— 1 



where cjj,c„ are the spinless fermion (SF) creation and 
annihilation operators, respectively. In fact, the Hamil- 
tonian (2) is transformed to a free spinless fermion model 

as 



(4) 



U = -\ y )Z ( C « C ™+! + C n+l C « 
n 

J* , V 

+ (4 C «+2 + 4+2 C «) 



We can diagonalize the Hamiltonian (4) by means of 
the Fourier transformation to obtain 



where 



H = ^E(k)c\k)c{k) 



E(k) = -J (cos(fc) - ^ cos 2kj . 



(5) 



(6) 



and a = J* /J. 

In what follows we assume that J > and a > 0. 
However, since the unitary transformation 



QZ QZ 



(7) 



changes the sign of the nearest neighbor exchange con- 
stant J — > — J, one can easily reconstruct the features 
of the system for J < using the transformation (7). 
Moreover, since the sign of the three spin coupling term 
is changed by the time reversal transformation 



qy _> qy qx,z 



(8) 



the properties of the system for a < can be eas- 
ily obtained from the ground state phase diagram of 
the system at a > using the transformation (8). 
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FIG. 2. The spectra for different values of the parameter 
a. For given a the corresponding curve is scaled so as to give 
the same bandwidth equal to 2. 

In the thermodynamic limit (N — -> oo), the ground 
state of the system corresponds to the configuration 
where all the states with E{k) < are filled and E{k) > 
are empty. For a < 0.5, the band minimum is at 
k = (see Fig. 2). If a > 0.5, the band minima 
±k min move away from k = and satisfy the relation 
cos(fc m i„) = l/2a, while a band maximum at k = is 
£■(0) = J (a/2 — 1). Therefore, for a < 2 there are only 
two Fermi points at 



arccos 



(l - Vl + 2a 2 ) /2a 



(9) 



and the ground state corresponds to the configuration 
when all states with |fc| < kp are filled. At a = a c = 2 
the band maximum at k = reaches the Fermi level (see 
Fig. 2) and at a > a c E(0) > and therefore, two 
additional Fermi points appear at 



±kp~ = arccos 



(l + \J\ + 2a 2 ) /2a 



(10) 



In this case the ground state corresponds to the con- 
figuration when all states with kp < \k\ < kp are filled. 

As we show in this paper, the change of the topol- 
ogy of the Fermi surface of the equivalent SF model 
for a > a c = 2 corresponds to the phase transition in 
the ground state of the spin system. This transition is 
the second order quantum phase transition, which at the 
transition point is characterized by 

• non-monotonic behavior of the ground state energy 
Eq (a) as a function of the parameter a . In partic- 
ular, the first derivative of the ground state energy 
with respect to the parameter a shows a kink, while 
the second derivative is divergent at a = a c ; 

• magnetization of the system as a function of the 
parameter a shows the non-monotonic and nonan- 
alytical behavior at a = a c ; 
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• the critical index characterizing the power-law de- 
cay of the transverse spin-spin correlation function 
is changed at a — a c ; 

• At low temperature (T <C J) the heat capacity, 
magnetization and magnetic susceptibility of the 
system show a well pronounced anomalous behav- 
ior at a ~ a c ; 

This quantum phase transition, despite its obviously 
rich nature, can be studied using the exact solution of 
the model (2), given in terms of Fermi gas of spinless 
particles (5). 



III. PHASE TRANSITION IN THE GROUND 
STATE 

The ground state energy of the system is given by 

£o(a) = %-\ £(k)dk (11) 

27r J {A} 

where the integration region {A} = [—hp, hp] for a < 2 



kp, k F ] 



U [kj.,kp] for a > 2. 



and {A} = [- 

In Fig. 3 we plotted the ground-state energy of the sys- 
tem as a function of the parameter a. At a = a c we 
observe (see inset in Fig. 3) the singularity in the behav- 
ior of the generalized stiffness 

f(a) = -d 2 £ a (a)/d 2 a. 

The singularity in £(cv) we attribute to the 
second-order phase transition in the ground- 
state of the system with the increase of the 
three-spin coupling constant J* = aJ. 




"0 0.5 



FIG. 3. The ground state energy as a function of the pa- 
rameter q. The inset shows the generalized stiffness £ as a 
function of the parameter a 



A. Magnetization 

Most clearly, the change of the properties of the sys- 
tem at the transition point is seen in the non-monotonic 



behavior of magnetization. At T — the magnetization 
(per site) of the system is given by the number of SF in 
the ground state 



™'(«) = lM = 
kp/ir- 1/2 at (a < 2) 

(kp - fc£)/7r- 1/2 at (a>2) 



(12) 



Using Eqs (9)-(10) we found that the ground-state of 
the system is singlet with m z — only at a = (the XY 
model) and in the limiting case a = oo. For arbitrary 
nonzero a, the magnetization of the system is finite. It 
is a monotonically increasing function of the parameter 
a in the range of < a < 2 (see Fig. 4). It reaches its 
cusp-type maximum at a = 2 with m z = 1/6 and then 
monotonically decreases to zero with a — > oo. 

The derivative of magnetization with respect to the 
parameter a (see Inset in the Fig. 4) is 



dm z (a)/da 



constant 
l 

\/ol — a c 



at (a — a c — > 0-) 
at (a — a c — > 0+) 




FIG. 4. The ground state magnetization as a function of 
the parameter a. The inset shows the derivative of magneti- 
zation with respect to the parameter a. 



B. Correlation functions 

To clarify the symmetry properties of various phases 
realized in the ground state of the system, we calcu- 
late the longitudinal and transverse spin-spin correlation 
functions 



JC*(r) = (S*S* +r ), 
1C tr (r) = (SZSZ +r ) = (S«S y n+r ). 



(13) 
(14) 



The spin-spin correlations in the case of ordinary XY 
chain (a = 0) were studied in the classical paper by Lieb, 
Schulz and Mattice 6 . The longitudinal spin-spin corre- 
lation function is 



lC*(r) 
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and the transverse correlation function decays more 
slowly 



JC tr {r) 



-1/2 



For further results on these correlation functions, see the 
Rcf. 11 . 

In our calculations of correlation functions we follow 
the route used in Ref. . In this approach the correlation 
functions are given in terms of Tocplitz determinants. 
Using the Jordan- Wigner transformation (3) one can eas- 
ily found that 



IC z (r) = del 



3(0) g{r) 
9(r) 9(0) 

H5 2 (0)-.9 2 M). 



(15) 



and 



K tr (r) = detG = 
5 (1) 3 (2) ... g(r) 

5(0) g(l) 



5 (2 -r) 5 (3 -r) 



g(r - 1) 



(16) 



Here 



g(r) = — sin(/c F r) - S 0r 



at a < a Cl and 



g(r) = — [sin(kpr) - sin(fc^r)] - S 0r 



(17) 



(18) 



at a > a c . 

After straightforward calculations we obtain that for 

a < a c 



K, z {r)=m 



1 sin 2 (fc F • r) 
7T 2 r 2 



(19) 



and for a > a c 



^ „ 1 (sin(fc p • r) — sin(fc± • r)) , 

As we observe, the quadratic decay of the longitudi- 
nal correlations remains unchanged at the transition 
point. However, the additional oscillations in the lon- 
gitudinal correlation function ~ cos(2/c F • r) and ~ 
cos {{kp ± hp) ■ r) associated with the presence of two 
Fermi points are clearly seen. 

The transverse correlation undergoes more dramatic 
change. We obtain that, for a < a c 



rtr M _ A. , B ■ C0S t 2fc F • r) 

' V V) r l/2 ^ r 5/2 



(21) 



and for a > a c . 



^ = Bi cos(fc F r + (^_)) + Ci cos (k+r + ip + ) 



+ 



+ 



r r 
D\ cos ((2kp — kp)r + ipi) 

Ei cos ((2kp — kp)r + ip 2 ) 



(22) 



Here A = A(a), ...E\(a) are the smooth functions of the 
parameter a. 

Below the gapless phase with the power-law decay of 
spin-spin correlations given by Eqs. (19) and (21) is 
called the Spin Liquid I (SL-I) phase, while the gapless 
phase with the power-law decay of spin-spin correlations 
given by Eqs. (20) and (22) - the Spin Liquid II (SL-II) 
phase. 

The behavior of the correlation functions given by Eqs. 
(19)-(22) perfectly fits into the usual conformal theory re- 
sults for the systems with several gapless excitations 12 ' 13 . 



C. Emptiness formation probability 

The very important quantity which characterizes a 
quantum spin system in the spin-liquid phase is the so 
called emptiness formation probability (EFP) 14 , i.e. the 
probability to find a ferromagnetic string of the length 
"n" in the spin liquid ground state 



n 

P(n) = (GS\l[(S* + -)\GS) 



(23) 



To calculate the EFP we follow the route developed in 
the recent paper by M. Siroishi, M. Takahashi and I. 
Nishiyama 15 . In the SF representation the EFP is de- 
scribed in terms of the following determinant 



P(n) = 



fll fl2 ■ ■ ■ fin 
/21 /22 • • • f2n 



fnl fn2 • • • fn 



(24) 



where 



f nm = (GS\cl c m \GS) 
= l/(2ir) j A dq e~xp(—iq(n — m)) . 



(25) 



For a < a c 



ir(n — to) 



and for a > a c 
1 



ir(n — to) 
2 



sin(fc F (n — m) /2) cos(fc F (n — to) /2) 

(26) 



[sin(fc F (n — to)) — sin(fc F (n — to))] 



sin([fc F — kp](n — m)/2) 



ir(n — to) 
x cos([fc F + kp](n - to)/2) 



(27) 
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FIG. 5. The emptiness formation probability P(n) at n = 2 
and n — 8 as a function of the parameter a. 

Using Eqs. (24)- (27) we calculated the exact values of 
EFP P(n) for n = 1, 20. The results of these calcula- 
tions are presented in the Table 1 (see Appendix A). In 
Fig. 5 we plotted the EFP P(n) as a function of the pa- 
rameter a for n — 2 and n = 8. As it is clearly seen from 
this figure the EFP shows the nonmonotonic behavior at 
a = a c . 

To consider whether QPT leads to the change in the 
asymptotic (n — > oo) behavior of EFP we derived the 
analytical expressions for EFP at a < a c and a > a c 
which fit the numerical data. 

At a < a c , (in full agreement with the results obtained 
in the Ref. 15 ) we obtain the following expression for EFP 



P<(n)=C(cos^ 



-1/4 



(28) 



where C = 0.6450. 

At a » a c we obtain the following expression for EFP, 
which fits our exact data 



P»(n) = 



sin 



sin 



k+ 



A 1 + (-l) n B 1 



n 



1/2 



(29) 



where A\ = 0.659 and B\ = 0.054. In the intermediate 
regime, for a > a c the fitting formula is more compli- 
cated and is presented in Appendix B. As we see, the 
QPT in the system manifestly changes the EFP. 



D. Order parameter 

The QPT at a = a c from the SL-I to the SL-II phase 
requires an Order Parameter for complete description. 
As the order parameter we introduce the following quan- 
tity 



77(a) = C E xtxY - JCsl-i 



where 



C 



ExtXY = 



EP(n) 



(30) 



(31) 



is the average length of the ferromagnetic string in 
the ground state of the extended XY chain with cou- 
pling a and, respectively, the magnetization m z (a), 
while Csl-i is the same quantity calculated in the SL- 
I phases with the same value of the magnetization. 




FIG. 6. Order Parameter as a function of the parameter a 

Straightforward calculations give (see Fig. 6 ) 
a?; (a — a c — > 0_) 



77(a) 



II 

y/a — a c at (a — a c ) — > 0+) 



(32) 



Thus, the order parameter exhibits the standard mean- 
field type behavior. 



IV. THERMODYNAMICS 

The QPT with respect to the parameter a in the 
ground-state of the system (2) is smeared out by thermal 
fluctuations at finite temperatures. Thus, the thermo- 
dynamic quantities have no singularities but, in the low- 
temperature region T <C J they show a well pronounced 
anomalous behavior at |a — a c | <C T/J. 

To investigate the thermodynamics of the system let 
us first calculate the free energy. In the thermodynamic 
limit the free energy of the system is given by 



F = -T j dujp(Lu) In [2 cosh (w/2T)] . 



(33) 



The SF density of states p(uj) is nonzero and equal to 



p(uj) = p (uj) 



(34) 



for a < 2 and for a > 2 and E(0) < uj < W max . For 
a > 2 and W min < w < E(0) 
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p(w) = p + p + (uj). 



(35) 



Here W min = min{-J(2 - a)/2,-J{2a 2 + l)/4a}, 
W maK = J(a + 2)/2, and 



1 2a 



:, (36) 



where 



G a (oj/J) = y/T+ 2a 2 + 4a(ui/ J) . 



(37) 



Knowing the free energy it is easy to obtain the ex- 
pression for the entropy 



,S = /_ 00 co ^H[ln(2cosh(^)) 
"(#) tanh(^) 



(38) 



Below in this Section we study the anomalies in the 
low-temperature behavior of different thermodynamic 
quantities, caused by the presence of the QPT in the 
ground state at a ~ a c . 



A. Specific heat 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
T/J 

FIG. 7. The specific heat of the system as a function of the 
parameter T/J for different values of the parameter a. The 
inset shows the specific heat of the system at T/J <C 1. 

In Fig. 7 we plotted the heat capacity of the extended 
XY model for different values of the parameter a. Since 
the bandwidth of the system depends on the parame- 
ter a, in Fig. 7 we plotted the results for the models 
with normalized bandwidth Wr = 23 . In the inset we 
present the low-temperature [T <C J) heat capacity. The 
anomaly in the behavior of heat capacity at a = a c = 2 
is clearly seen. 



The anomalous temperature dependence in the vicinity 
of the QPT is seen most clearly in the low-temperature 
behavior of the specific heat. The heat capacity of the 
system is given by 



C 



J —oo 



cosli (w/2T) 



(39) 



Away from the critical point, the density of states in 
the vicinity of the Fermi level p~(0) at a <C a c and 
p~ (0) + p + (0) at a » « c , respectively, is the constant 
of the order of <~ 1/J. Therefore, in these cases the spe- 
cific heat of the system 

C(T)^ 7 (T/J) 

with 7 <~ 1. 

However, in the vicinity of the critical point, at 
\ot — a c \ < \co/J\, 



p+{u) ~ 



1 



7ry 6 J\l 



(40) 



and therefore, in this case the heat capacity of the sys- 
tem exhibits the anomalous square root dependence on 
the temperature 



C = 'n-y/f/J + 'Y-(T/J). 



(41) 



B. Magnetization and magnetic susceptibility at 

T + 0. 

At finite temperature the magnetization of the system 
is proportional to the average number of spinless fermions 

1 [°° 

mZ ak T ) = -T > \ duj p(uj) tanh (u /2T) , (42) 

2 J-00 



while the magnetic susceptibility is given by 
1 f°° 1 

x(T) = i ^y_ oo ^M cosh2(a;/2T) . 



(43) 



Let us first consider the limiting case of small temper- 
atures T/J C 1. In this case it can be easily obtained, 
that for a <C a c 



m' a (T) = m* (0) (^1 - - ■ (T/jy + 0((T/Jf) j , (44) 

Xa{T) = Xa {0) (l + y 2 ■ (T/J) 2 + G((T/J) 4 )) , (45) 
and for a 3> ot c 

m* a (T) = m* (0) (1 + 2tt 2 • (T/J) 2 + 0((T/jf)) , (46) 

Xa(T) = xa(0) (l + ^ • (T/J) 2 + 0((T/Jf)\ . (47) 



6 



• The low temperature (T/J <C 1) expansion shows 
that magnetization of the system decays with tem- 
perature at a < a c , and magnetization of the sys- 
tem increase with the increase of the temperature 
at a > a c . 

• The magnetic susceptibility of the system is finite 
both for a <C a c and a » a c and shows the in- 
verse square root dependence on temperature for 
\a — a c \ -C T/J. 

In the case of coupling close to the critical value, at 
\a - a c \ < T/J 

<(T) = m* (0) [l - /3 1/2 • (T/jf 2 - (3 3/2 ■ (T/J) 3 / 2 

- 02 ■ (T/J) 2 + 0((T/ J) 5/2 )] , (48) 
X(T) = /?_i/ 2 • (TJ)- 1 ' 2 + p-(u = G) + 0(T/J) 2 . (49) 
Here 

V2- 1 



: C(l/2) >0, 

V 1Z7T 

:! ^C(3/2)>0, 



96V6tt 

7T 



Pl/2 
/?3/2 

ft 

0-1/2 

and ((z) is the Riemann's zeta function. 



81 V3 ' 

-^=(2>/2-l)C(-l/2)>0, 



(50) 

(51) 
(52) 

(53) 




FIG. 8. Magnetization of the system as a function of tem- 
perature for different values of the parameter a. The inset 
shows the magnetic susceptibility of the system as a function 
of temperature for different values of the parameter a. 

Let us now consider in detail the problem of anomalous 
temperature dependence of magnetization for a > a c . In 
Fig. 8 we plotted the magnetization of the system as a 
function of temperature for three different values of the 
parameter a = 1, 2 and 4. From Fig. 8 it is clearly seen 
that at a = 4 the magnetization reaches its maximum at 
T = 0.15 J and then smoothly decays with the increase 
of temperature. 



To gain more insight into the behavior of magnetiza- 
tion in Fig. 9 we plotted the magnetization as a function 
of the parameter a for different values of the temperature. 
As it is clearly seen from Fig. 9, M(T) < M = M(T = 
0) at arbitrary T when a < a c . But at a > a c there does 
exist the crossover temperature T c , and so M(T) > Mq 
for < T < T c and M(T) < M for T > T c . We cal- 
culated T c as a function of the parameter a (see inset 
in Fig. 9) and found that T c monotonically increases 
with the increase of a from its minimum value T c = at 
2 up to T c = 0.186 J at a 



a = a c 
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FIG. 9. Magnetization of the system as a function of the 
parameter a at different values of temperature. 

The source of such unusual behavior is the struc- 
ture of the excitation spectrum: at a =/= the ground 
state of the system is characterized by the magnetization 
M^. At a < a c , the excited states with magnetization 
M = M^±m have equal energy for arbitrary m = 1, 2, .... 
Therefore, for a < a c at arbitrary finite temperature 
M a {T) < Mq. On the other hand, for a > a c the excita- 
tion spectrum of the system exhibits a clear anisotropy: 
at least some lowest excited states with magnetization 
M = Mq + m (with m = 1, 2..., m ) are far below their 
counterparts with magnetization M = Mq — m. There- 
fore, at sufficiently low temperatures, where the states 
with magnetization M > are occupied mainly, the 
total magnetization of the system increases. However, 
with further increase of the temperature, the difference 
in occupation of states with M = M° ± m reduces, and 
finally, at T > T c , the magnetization of the system drops 
below its ground state value. 



V. PROPERTIES OF THE SYSTEM IN THE 
PRESENCE OF EXTERNAL MAGNETIC FIELD 

In this Section we consider the properties of the ex- 
tended XY spin chain (2) in the presence of an external 
magnetic field. We consider the general case, where the 
field dependent part of the Hamiltonian is given by 



(54) 
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In terms of spinless Fcrmions the Hamiltonian of the 
extended XY model in the presence of an external mag- 
netic field reads 

J^/ + t 



w = -g 1] (4 C «+ 

n 

+ — XI (4 C "+2 + 4+ 2 Cn) 
n 

-£(#0 + (-im) (4^-^ 



(55) 



Below in this Section we consider separately the ground 
state phase diagram of the extended XY spin chain in 
the case of: 

• applied uniform magnetic field H$ = 0; 

• applied staggered magnetic field H = 0; 

• applied magnetic field acting only on the spins on 
the even sites H = H$ = H. 

A. Uniform magnetic field 

In the case of applied uniform magnetic field (Hs = 
0, Hq = hJ) the diagonalization of the Hamiltonian is 
straightforward and gives 



H = ^E 1 (k)c\k)c(k) 



(56) 



where 



J E 1 (fc) = -j(/i + cos(fc)-|cos2fc) . (57) 

The phase diagram of the model in the case of applied 
uniform magnetic field is presented in Fig. 10: 

• For hJ < W m in and hJ > W max the system is in 
the ferromagnetic phase . 

|a) the system is in 



For - (1 
the SL-I phase. 



\a) < h < (l 



• For a > 0.5 and - (2a 2 + l) /4a < h< - (l - \a) 
the system is in the SL-II phase. 




D 0.5 



FIG. 10. The ground state phase diagram of the model in 
the case of applied uniform magnetic field. 



B. Staggered magnetic field 

In the case of applied staggered magnetic field (Hq = 
Jh, Hs = Jhg) it is useful to introduce two types of 
spinless fcrmions determined on the odd and even sub- 
lattices, respectively 



Cln-\ — a n -i/2 



and 



C2n = b n . 



(58) 



In terms of these new particles, the Hamiltonian in the 
momentum space reads 



U = X \^{e a {k)a\a k + e b (k)b 

k 

£ab(k) (a{b k + b\a^j 



where 



+ 



e a (k) = — cos(2fc) - h Q + h s , 



(59) 



e b (k) 



■ cos(2fc) — ho — hs , 



(60) 



- Jcosfc . 

Using the standard Bogoliubov rotation 



"/ = <*# I ^0fc J a k + sin (^0 k ) I,, 



b k 



sm | -0 k ) a k + cos \-0k ) Pk 



where 



tan(0 fe ) = 



2s ab (k) 



e a {k) - e b {k) 
we finally obtain 

7i = J2 {E-{k)a\a k + E + {k)0[f3k) , 
k 

where 

E± (k) = J cos 2k - h ± ^h 2 s + cos 2 kj . 



(61) 



(62) 



(63) 



(64) 



Let us first consider the case of the staggered magnetic 
field (Ho = 0). The ground state phase diagram of the 
model in this case consists of the following four sectors 
(see Fig. 11 ): 

• At \hs\ > \ct the system is in the long-range or- 
dered Neel anti-ferromagnetic phase . 
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There is a gap in the spin excitation spectrum. The lon- 
gitudinal spin-spin correlation function is 

JC z (r) = {S z n S z n+r ) ~ {-Iff {a) 

where f{a) is the constant of the order of unity for given 
a. 



• At a < 2 and \hs\ < \a as well as for a > 2 and 



q,2_ 4 



< \hg\ < ha the system is in the 5X-/phasc 




At a > 2 and |ft 5 | < 
SL-II phase. 



Q 2 -4 



the system is in the 




FIG. 11. The ground state phase diagram of the model in 
the case of applied staggered magnetic field. 



C. Mixed magnetic field 



In this subsection we consider a rather special case, 
where the mixed magnetic field with equal strength of 
the uniform and staggered components Hq = Hs = hJ. 
This corresponds to the case, where the magnetic field of 
the 2hJ strength is applied to the spins on the even sites, 
while the spins on the odd sites experience no magnetic 
field. 

The spectrum of the system in this case is given by 



Ef (k) = J (| cos 2k - h ± ^Jh 2 + cos 2 fc) . (65) 



FIG. 12. Phase diagram for Uniform Magnetic field + al- 
ternating Magnetic field 

The ground state phase diagram of the model (59) in 
the case of applied mixed magnetic field consists of the 
following four sectors (see Fig. 12) 

• At a < \/2 and — 4 ~° < h < —\ot the system 
is in the long-range ordered Need antiferromagnetic 
phase ; 



At a > \[2 and — \a. < h < — 4 4 ° 2 the system 
in the SL-II phase; 



is 



At h > max{— \a; — 4 4 ° 2 } the system is in the 
SL-I phase; 



At h < min{— \a; — 4 4 " 2 } the system is in the SL- 
III phase; 



The new SL-III phase corresponds to the gapless phase 
characterized by qualitatively different behavior of the 
long-range spin-spin correlations for the spins located on 
the even and odd sites. The calculation of the longitudi- 
nal spin-spin correlations between the spins on the even 
sites gives 



£o, 2m (r = 2m) = (™r") 2 , 
between the spins on the odd sites 

^, 2m+1 (r = 2m) = (mf) 2 -^ 
and between the spins on even and odd sites 



q 2m (r = 2m-l)=m 



even ^odd 



(66) 



(67) 



(68) 



respectively. Here mf en < and m° dd > is the mag- 
netization per site of the even and odd sub-lattices, re- 
spectively, which is calculated at the given value of the 
coupling a. The parameter B is a smooth function of the 
coupling constant a and at the given a it is the constant 
of the order of unity. For example, at h = —0.4 and 
a = 1.5 mf en = -0.3406 and mf d = 0.1690. There- 
fore, in the SL-III phase, we have antiferromagnetically 
coupled two ferromagnetically ordered sub-lattices. 
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In the SL-III phase the transverse correlations also 
show the unusual behavior. In particular, the transverse 
spin correlation between the spins on the even sites is 

/C*y 2m (r = 2m) = 0, (69) 

between the spin on the odd sites 

(_l)m 

Ki,2 m+ i(r = 2m) * (70) 
and between the spins on the even and odd sites 

JCZ 2m {r = 2m-l)c y -J r . (71) 



VI. CONCLUSION 

In conclusion, we investigated the ground state phase 
diagram and the low-temperature thermodynamics of the 
extended XY model with nearest-neighbor exchange (J) 
and three-spin interaction J* = aJ. At J* ^ the 
spin system is characterized by the broken time rever- 
sal symmetry and for arbitrary a ^ 0,oo shows a fi- 
nite magnetization in the ground state M° (a) . We have 
shown that with the increase of the three-spin coupling 
the Quantum Phase Transition from the Spin-Liquid-I 
into the Spin-Liquid-II phase takes place. In the Spin- 
Liquid-I phase, at a < a* c — 2 the properties of the 
extended XY chain are similar to those of the standard 
XY model in the presence of such an "effective" mag- 
netic field that ensures the same value of the GS mag- 
netization M . However, in the Spin-Liquid-II phase, at 
a > a c , the behavior of the magnetization, magnetic sus- 
ceptibility, emptiness formation probability and different 
spin-spin correlation functions is qualitatively different 
and could not be considered by the simple effect of an 
"effective" magnetic field. The anomalous behavior of 
the heat capacity of the system, entropy, magnetization 
and magnetic susceptibility connected with the presence 
of QPT in the ground state has been shown. The ground 
state phase diagram of the model in the case of applied 
uniform and/or staggered magnetic field has been also 
obtained. 
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VIII. APPENDIX A 



In this appendix we present the exact values of the EFP P(n) for different values of the parameter a, obtained by 
evoluation of the Toeplitz determinant (24). 



TABLE I. The emptiness formation probability P(n). 







a = 


: 




a = 


: 1 




a = 


2 


a = 


4 




a = oo 






(m z ) 


= 


(m z ) = 


0.1193 


(m z ) = 


= 1/6 


(m z ) = 0.0213387 


(m z ) = 


n 


= 1 


0.5 




6.1928 


•10 


-1 


6.6667 


•io- 1 


5.2134 


■10 


-1 


0.5 


n 


= 2 


1.4868 


■10" 


-1 


2.9576 


•10 


-1 


3.6845 


•io- 1 


2.6455 


■10" 


-1 


0.25 


n 


= 3 


2.3679 


■10" 


-2 


1.0251 


•10 


-1 


1.6136 


•io- 1 


8.4078 


■10 


-2 


7.4339 -IO" 2 


n 


= 4 


1.9453 


■10" 


-3 


2.5041 


•10 


-2 


5.4640 


•io- 2 


2.5783 


■10 


-2 


2.2105 TO- 2 


n 


= 5 


8.1263 


■10" 


-5 


4.2525 


•10 


-3 


1.4127 


•io- 2 


4.5135 


■10 


-3 


3.5205 TO- 3 


n 


= 6 


1.7152 


■10" 


-6 


4.9869 


•10 


-4 


2.7702 


•10~ 3 


7.6228 


■10 


-4 


5.6069 TO- 4 


n 


= 7 


1.8232 


■10" 


-8 


4.0241 


•10 


-5 


4.1052 


•io- 4 


7.0945 


■10 


-5 


4.6062 TO- 5 


n 


= 8 


9.7402 


io- 


11 


2.2298 


•10 


-6 


4.5877 


■IO" 5 


6.3746 


■10 


-6 


3.7841 TO- 6 


n 


= 9 


2.6122 


io- 


13 


8.4729 


•10 


-8 


3.8609 


•10~ 6 


3.1110 


■10 


-7 


1.5808 TO- 7 


n 


= 10 


3.5137 


io- 


16 


2.2060 


•10 


-9 


2.4448 


•io- 7 


1.4664 


■10 


-8 


6.6036 TO- 9 


n 


= 11 


2.3691 


10" 


19 


3.9329 


-LU 


-11 


1.1641 


•io- 8 


3 79Q1 




-10 


1.3938 TO- 10 


n 


= 12 


8.0036 


10" 


23 


4.7992 


10" 


-13 


4.1656 


io- 10 


9.1615 


10" 


-12 


2.9420 TO- 12 


n 


= 13 


1.3543 


io- 


26 


4.0070 


10" 


-15 


1.1200 


io- 11 


1.2099 


10" 


-13 


3.1272 TO- 14 


n 


= 14 


1.1475 


10" 


30 


2.2884 


10" 


-17 


2.2619 


io- 13 


1.5439 


10" 


-15 


3.3240 TO- 16 


n 


= 15 


4.8677 


10" 


35 


8.9382 


10" 


-20 


3.4304 


io- 15 


1.0568 


10" 


-17 


1.7758 TO- 18 


n - 


= 16, 


1.0336 


10" 


39 


2.3871 


10" 


-22 


3.9063 


io- 17 


6.9905 


10" 


-20 


9.4872 TO- 21 


n 


= 17 


1.0984 


io- 


44 


4.3587 


10" 


-25 


3.3395 


io- 19 


2.4771 


10" 


-22 


2.5443 TO- 23 


n = 


= 18 


5.8409 


10" 


50 


5.4406 


10" 


-28 


2.1431 


io- 21 


8.4825 


■10" 


-25 


6.8235 TO- 26 


n 


= 19 


1.5534 


10" 


55 


4.6420 


10" 


-31 


1.0323 


1Q -23 


1.5538 


■10" 


-27 


9.1783 TO- 29 


n 


= 20 


2.0604 • 


io- 


61 


2.7070 • 


10" 


34 


3.7316 


1Q -26 


2.7538 


■10" 


-30 


1.2346 TO- 31 
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IX. APPENDIX B 

In this appendix we present the expression for the EFP, which fits our exact data for a > a c : 

Df , [. f kp-k+ \ . f k F +k+ \]^ f A 1 + (-l) n B 1 Ci+ D lC os(k+n) \ 

For different values of the parameters a the fitting values of the coefficients Ai, ...D\ are given in Table II. 



TABLE II. The coefficients of the fitting formula (72) for the emptiness formation probability P(n) for different values of 
the parameter a > a c . 



a 


A 1 








2.005 


0.155 


0.009 


0.302 


0.314 


2.01 


0.204 


0.011 


0.370 


0.194 


2.05 


0.331 


0.018 


0.350 


0.072 


2.1 


0.351 


0.021 


0.322 


0.063 


2.2 


0.572 


0.032 


0.150 


0.021 


2.3 


0.673 


0.037 


0.064 


-0.003 


2.4 


0.683 


0.039 


0.045 


-0.006 


2.5 


0.685 


0.041 


0.036 


-0.006 


3 


0.682 


0.045 


0.012 


-0.005 


4 


0.675 


0.050 


-0.0002 


-0.004 


10 


0.662 


0.053 


-0.006 


-0.003 


oo 


0.661 


0.054 
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